Mon. Not. R. Astron. Soc. 000, [Tb?? () 



Printed 30 January 2012 



(MN IaT E X style file v2.2) 



Intermediate- Mass Black Holes in Globular Clusters 



o 

(N 



Yu-Qing Lou 1 ' 2 ' 3 and Yi-Hong Wu 1 

1 Department of Physics and Tsinghua Center for Astrophysics (THCA), Tsinghua University, Beijing 100084, China; 

2 Department of Astronomy and Astrophysics, the University of Chicago, 5640 South Ellis Avenue, Chicago, IL 60637, USA; 

3 National Astronomical Observatories, Chinese Academy of Sciences, A20, Datun Road, Beijing, 100021, China 



Accepted 2012 January 26. Received 2012 January 7; in original form 2011 September 



Of 

6 



ABSTRACT 

There have been reports of possible detections of intermediate-mass black holes 
(IMBHs) in globular clusters (GCs). Empirically, there exists a tight correlation be- 
tween the central supermassive black hole (SMBH) mass and the mean velocity dis- 
persion of elliptical galaxies, "pseudobulges" and classical bulges of spiral galaxies. 
We explore such a possible correlation for IMBHs in spherical GCs. In our model of 
self-similar general polytropic quasi-static dynamic evolution of GCs, a criterion of 
forming an IMBH is proposed. The key result is Mbh = Ca 1 ^ 1 ^ 11 ^ where Mbh is the 
IMBH mass, a is the GC mean stellar velocity, £ is a coefficient, and 2/3 < n < 1. 
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1 INTRODUCTION 

For star-forming molecular clouds, collapsed massive dense 
gas cores eventually lead to luminous new-born stars burn- 
ing nuclear fuels. Analogously but on larger scales, we specu- 
late that grossly spherical core-collapses of globular clusters 
(GCs) could also cause something singular around the dense 
centre; such central singularities may form IMBHs. Observa- 
tions of GC cores indicate that the central concentrations of 
nonluminous materials are likely due to IMBHs (Bahcall & 
Ostriker 1975). Mass accretions onto IMBHs were proposed 
to power GC ULX sources (e.g. Farrell et al. 2009). 

On much larger scales further, observations of galaxies 
reveal a strong correlation between the mass Mbh of the cen- 
tral SMBHs and the mean velocity dispersion a of the host 
stellar bulge in the form of log(MBH/M©) = e + S log(a/o"o) 
where e and 5 are two coefficients and an is a velocity disper- 
sion normalization (e.g. Tremaine et al. 2002). We naturally 
expect a similar Mbh — a power-law relation for GCs. 

Evidence for central IMBHs in GCs have been debated 
extensively and their possible existence bears important con- 
sequences for both the formation and evolution of GCs (e.g. 
Grindlay & Gursky 1976; Maccarone et al. 2007, 2008; Za- 
harijas 2008). The central brightness excesses observed in 
several GC cores (e.g. Diorgovski & King 1986) are com- 
patible with the presence of central IMBHs. Properties of 
IMBHs correlate with various properties of GCs, including 
stellar density profiles, central stellar dynamics, luminosities 
of GCs, the mass-to-light ratios, surface brightness, rotation 
amplitudes, position angles, dark matter densities and the 
mean stellar velocity dispersion a (Gebhardt et al. 2000; 
Zheng 2001; Ulvestad et al. 2007; Zepf et al. 2007; Bash et 



al. 2008; Noyola et al. 2008; Zaharijas 2008). These empirical 
correlations may shed light on the origin and evolution his- 
tories of IMBHs and their host GCs. Among such observed 
relations, Mbh and a correlate tightly (e.g. Gebhardt et al. 
2002; Safonova & Shastri 2010). 

Self-similar solutions for general polytropic hydrody- 
namics of a self-gravitating fluid with spherical symmetry 
were constructed recently. Asymptotic behaviours of novel 
quasi-static solutions in a single polytropic fluid have been 
revealed by Lou & Wang (2006, 2007) and was used to model 
rebound (MHD) shocks in SNe. Such solutions were applied 
to clusters of galaxies (Lou et al. 2008) for possible galaxy 
cluster winds. In this Letter, we invoke such quasi-static 
solutions to model dynamic evolution of host GCs and for- 
mation of IMBHs and to establish Mbh — a power laws. 

Various aspects of GCs have been studied extensively 
(e.g. Benacquista & Downing 2011 and extensive references 
to excellent reviews therein). We focus on the quasi-static 
self-similar GC dynamic evolution (say, induced by the 
gravothermal instability) in the late phase of the pre-collapse 
regime; such asymptotic solution for GC evolution leads to 
diverging mass density at the center. Physically, as the inner 
enclosed core mass becomes sufficiently high within a radius 
comparable to its Schwarzschild radius, an IMBH forms in- 
evitably (e.g. through mergers of stellar mass black holes 
or runaway collisions and coalescence or merging of stars to 
form supermassive stars and to trigger subsequent e* pair 
instabilities therein). After forming such a central IMBH in 
the core, pertinent post-collapse mechanisms continue to op- 
erate: e.g., mass segregation maintains more massive stars 
around the central IMBH, while the 'binary heating' (e.g. 
Hurley et al. 2007) from primordial stellar binaries survived 
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the IMBH formation tends to resist further core collapse or 
to drive post-collapse oscillations. The N-body GC simula- 
tions by Hurley et al. (2007) of up to 10 s stars and initial 5% 
binaries may eventually reach a GC core binary frequency 
as high as 40% at the end of the core-collapse phase. Shown 
in their figure 3, this core binary frequency actually fluc- 
tuates between ~ 10% to ~ 40%. We would expect that 
an IMBH forms at the GC centre rapidly during the core- 
collapse phase but do not know exactly when. Physically, 
this IMBH would engulf the central stars and binaries at 
the epoch of IMBH formation. Such an explosive event may 
give rise to a powerful gamma-ray burst and a shock wave 
surrounding the GC centre. The slower relaxation, evolution 
and accretion then persist on a much longer time scale. 

GCs are close to spherical; in their dynamic evolution, 
lumpiness observed could either result from merging distur- 
bances and tidal disruptions or provide source of fluctuations 
that may be classified into acoustic modes, gravity modes 
and vortical modes on larger scales (Lou & Lian 2011). These 
perturbations may be unstable to trigger gravothermal in- 
stability. Analogous to nuclear burnings in a star to resist 
stellar core collapse, the 'binary heating' from primordial 
stellar binaries may delay core collapse in GCs (e.g. Mey- 
lan & Heggie 1997). As the source of 'binary heating' is 
exhausted during a GC evolution, the inner core collapse is 
inevitable. This ultimately leads to the formation of IMBH 
which may accrete materials from immediate environs. 



2 GC MODEL FOR A M BH o POWER LAW 

We adopt the same hydrodynamic perspective of Lou & 
Jiang (2008; LJ hereafter) for the large-scale spherical GC 
dynamic evolution. GCs are smaller and less massive than 
typical galactic bulges. The random stellar velocity disper- 
sion and the 'binary heating' from primordial stellar binaries 
in the core provide an effective pressure against the GC self- 
gravity. This may justify our fluid formalism for GC cores. In 
contrast to LJ, we emphatically focus on the reported ten- 
tative candidates of IMBHs and the properties of the host 
GCs. Our main goal is to extend the theory of LJ to GCs 
and examine whether IMBHs and properties of GCs can be 
sensibly fitted with data. By data comparisons, our results 
appear encouraging for such a physical connection. Our pre- 
dictions for IMBHs and host GCs can be tested by further 
observations. Meanwhile, we also compare with SMBHs in 
their host galactic pseudobulges. This may hint at a general 
validity of such a self-similar dynamic evolution in spherical 
self-gravitating systems. 

In spherical polar coordinates (r, 9, <j>), the nonlinear 
general polytropic hydrodynamic partial differential equa- 
tions (PDEs) of spherical symmetry are PDEs (1) — (4) of 
LJ with the same notations. The Poisson equation is auto- 
matically satisfied. As the bulk flow of stellar fluid is slow 
in GCs, we invoke the quasi-static self-similar solutions of 
Lou & Wang (2006). We introduce the transformation in the 
dimensionless independent variable x, 
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with K and n being two scaling parameters; here, u, p, P, 
and M are radial velocity, mass density, pressure, and en- 
closed mass respectively, while v(x), a(x), fi(x), and m(x) 
are respectively dimensionless reduced speed, mass density, 
pressure, and enclosed mass of x only. 

Substituting self-similar transformation JTj into non- 
linear PDEs (1) — (4) of LJ and defining q = 2{n + 7 — 
2)/(3n — 2), we derive two coupled nonlinear ordinary dif- 
ferential equations (ODEs) for a' and v' 

D(x,a,v)a = N\(x,a,v), D(x,a, v)v' = N2(x, a, v) , (2) 

where three functionals D, Ni and N2 are defined explicitly 
in Hu & Lou (2009) with zero magnetic field. An exact global 
static solution of eq ((2]) in physical dimensions, known as the 
singular general polytropic sphere, has u = 0, 
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where coefficient A 
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(3n-2)G 

, ■, -l/(n-3n«/2) 



(3) 



, , . This Mil 1 1 

2(2 — n)(3n — 2) 

tion serves as an asymptotic 'quasi-static' solution for small 
x; i.e. they are leading terms of v(x) and a(x) and there exist 
higher order terms for a self-similar asymptotic evolution. 

For such quasi-static self-similar hydrodynamic asymp- 
totic solutions at small x, we consistently presume 
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in two coupled nonlinear ODEs ((2)1 , and derive two nonlinear 
algebraic equations for the coefficients J, S and L, 
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where W = n 2 q/[2(2 — n)(3n — 2)]. Once the proper roots of 
S are known, coefficients J and L are related by eq ((4}; only 
one is free to choose. The existence of SPS solution (|3} and 
the requirement of > 1 constrain the parameter regime 
of such quasi-static solution. Oscillatory behaviours can also 
emerge (Lou & Wang 2006), which may be relevant to the 
interesting controversy of gravothermal and post-collapse os- 
cillations in GCs. For a sufficiently small W ^ 0, we can have 
two roots S > 1 from quadratic eq ([5]). It also happens for 
one root S > 1 and the other root S < 1. 

As a physical requirement for sensible similarity solu- 
tions of a general polytropic flow, both v(x) and a(x) ap- 
proach zero at large x. Thus for either x — > + or x -4 +00, 
the reduced velocity v — > 0, which means at time t, for either 
r — > + or r — > +00 the flow speed u —¥ 0, or at a radius r, 
when t is either short or long enough, the radial flow speed 
u — ¥ 0. This model describes a self-similar dynamic GC evo- 
lution towards a quasi-static configuration a long time later 
and may grossly fit relaxed spherical GCs. 

From the general polytropic EoS (LJ), the pressure is 

P = (47r) 7 - 1 G' J+7 - 1 (3n - 2) q K 1 - 3 ^ 2 p 1 M q 

and ol(t, t) = (7-P/p) 1//2 is the local stellar velocity disper- 
sion in a GC. Asymptotically as t — > +00, crt(r, t) becomes 



o"i(r) = 7 ' K 



1/2 ^l/(2n) o/2 ,( 9+7 -l)/2 («-l)/n 



(6) 
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To check against data, we derive the spatial average of ve- 
locity dispersion a in a GC. The GC boundary is taken as 
either radius r c where mass density p c is indistinguishable 
from the surrounding or the tidal radius. Within r c of a GC, 
the spatial average of stellar velocity dispersion o"l(V) is 

a — — 5— / <TL(r)47rr 2 dr = QK 1 ^ 2 
47rrg J 

= [3n 1+q/2 7 1/2 /(4n - l)](4nGp c ) (1 ~" )/2 A :inq/4 K 1/2 . 

We invoke a heuristic criterion of forming an IMBH in a GC 
(LJ). An IMBH mass Mbh is given by M BH = r s c 2 /(2G) 
where r s is its Schwarzschild radius and c is the speed of 
light. By solution (|3} and when 

(3^2)G r =rc/{2G), 

an IMBH forms with the Schwarzschild radius 

r = r s = [(3n - 2)c 2 /{2nAK 1/n )] n/{2n - 2) . 

Only those asymptotic quasi-static GCs with n < 1 can thus 
form central IMBHs (see fig. 1 of LJ). Consequently, 

M BH = [c 2 /(2G)][(3n - 2)c 2 /(2rM)] n/(2n " 2) K 1/(2 - 2n) , (7) 

or equivalently, the explicitly Mbh — <J power law 

n/(2-2n) / \ l/(l-n) 

= C<J 1/(1 - n) , 



10 J 



A/BH = ^ 



2nA 



(3n - 2)c 2 



where the exponent 1/(1 — n) > 3 since 2/3 < n < 1 (LJ). 

To validate our quasi-static model for the nine GCs with 
available observational data and references summarized in 
Tables 1 and 3, we fit a Mbh — a power law in Fig. 1. 
Applying the least-square criterion to the data of Table 1 
(e.g. Meylan & Mayor 1991; Safonova & Shastri 2010), we 
obtain Mbh = 4.1020 x 10 7 M o (cr/200 km s" 1 ) 3 ' 6251 with 
parameters {n, 7, p c } being (0.7241, 1.99, 4.40M e pc" 3 ). 
This p c value appears somewhat larger but grossly consis- 
tent with the data in rough orders of magnitude (e.g. Mey- 
lan 1987). Specifically, p c ~ O.87M pc" 3 for GC 47Tuc 
(Meylan 1988), p c ~ 0.091M o pc -3 for GC NGC6397 
(Meylan & Mayor 1991) and p c ~ O.42M pc" 3 for 
GC Gl (Meylan et al. 2001). No published p c values 
for NGC2808, M80, M62, NGC6388 and M15 are avail- 
able. The nine GCs 47Tuc, NGC2808, u Cen, M80, M62, 
NGC6388, NGC6397, M15, and Gl correspond to K = 
{1.2, 2.0, 7.7, 1.5, 2.1, 3.1, 0.22, 1.9, 5.8} x 10 19 cgs unit 
and r c = 13, 16, 32, 14, 17, 20, 5, 16, 28 pc, respec- 
tively in rough agreement with the estimated data in or- 
ders of magnitude. Harris (1996) reported r c of 47Tuc, 
NGC2808, u Cen, M62, NGC6388, NGC6397, and M15 to 
be 56, 43, 88, 18, 18, 11, and 64 pc, while Bahcall & Haus- 
man (1976) reported r c of M80 to be ~ 18 pc and Ma et al. 
(2007) estimated r c of Gl to be ~ 81 pc. 

We now consider a sample of SMBHs in galactic pseu- 
dobulges (e.g. Kormendy & Kennicutt 2004) summarized in 
Table 2. For the reasons in Hu (2008), we do not regard 
galaxy NGC3227 as containing a pseudobulge. In Table 2, 
cr e is defined by Gebhardt et al. (2000) and Tremaine et al. 
(2002) as the luminosity-weighted rms velocity dispersion 
within a slit aperture of length 2R e , where R e is the effective 
or half-light radius of a galactic bulge. Parameter as is the 
rms velocity dispersion within a circular aperture of radius 
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Figure 1. Mass of IMBH Mbh/A^q versus mean velocity 
dispersion cr/(200 km s _1 ) for the nine GCs (Tables 1 and 
3). The solid line is the least-square fit: Mbh = 4.1020 X 
lO 7 M (tr/2OO km s - 1 )3'625l witll parame ters {n, 7, p c } being 
(0.7241, 1.99, 4.40M Q pc" 3 ). 



Table 1. Nine GCs with their reported central IMBHs and mean 
stellar velocity dispersions a 

GC Name Other Name A/ B h(10 3 Mq) cr(km s" 1 ) 



NGC104 


47Tuc 


1 0+ ' 5 
l - u -0.5 


11.6 ± 0.8 


NGC2808 




2.7 


13.4 


NGC5139 


u> Cen 


30 ±4 


22.8 


NGC6093 


M80 


1.6 


12.4 ±2.5 


NGC6266 


M62 


3.0 


14.3 ± 0.4 


NGC6388 




r; 7+5.7 
°- ' -2.85 


18.9 ± 3.6 


NGC6397 




0.05 


4.5 ±0.6 


NGC7078 


M15 




14.1 ± 3.2 


Gl (M31) 


Mayal II 


18.0if : ° 


25.1 ± 1.7 



Re/8. Ferrarese & Merritt (2000) used "central stellar veloc- 
ity dispersion" a c . There is a relatiorQ erg = a c (8Rs, p / R e ) 04 
(J0rgensen et al. 1995), where i? ap ~ 2" (e.g. Davies et al. 
1987). Actually, the difference between er e and a c are much 
smaller than their errors (e.g. Hu 2008). 

A pseudobulge border is at radius r c where p 
reaches a value p c indistinguishable from the environs. 
With the least-square fit to Table 2 data, we obtain 
M BH = 2.4350 X 10 7 M Q (cr/200 km s -1 ) 3 7680 with pa- 
rameters {n, 7, p c } being (0.7346, 1.999, 1.8M@ pc" 3 ), 
respectively. This estimated p c appears grossly consis- 
tent with the data in the order of magnitudes, e.g. 
p c ~ O.IM© pc" 3 for the pseudobulge in the Milky Way 
(Lopez-Corredoira et al. 2005). No published p c values for 
other pseudobulge are available. The eight pseudobulges 
NGC1068, NGC2787, NGC3079, NGC3384, NGC3393, 
Circinus, IC2560, Milky Way correspond to values of K = 
{3.4, 5.8, 1.3, 3.5, 5.0, 0.85, 1.4, 2.6} x 10 21 cgs unit and 
Tc = 1.3, 1.7, 0.79, 1.3, 1.5, 0.64, 0.82, 1.1 kpc, respec- 
tively in rough agreement with the data in orders of mag- 



1 The prime distinguishes this approximation for as from the 
actual value of as and the ratio a' s / as may depend systematically 
on the velocity dispersion of a galaxy. 
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Table 2. A sample of SMBHs in pseudo-bulges of galaxies 



Table 3. Nine GCs with the reported total GC masses Mqq 



galaxy 


A/ BH (1O 8 M ) 


<r e (km s 1 ) 


(j c (km s - 1 ) 


NGC1068 a 


0.15 


165 ± 17 


165 ± 17 


NGC2787 6 


41+ 04 


210 


210 


NGC3079 

NGC3384 C 

NGC3393 


025+ ' 025 
u.uzo_ 013 

16+ 001 
u - 1D -0.02 

0.31 ±0.02 


146 ± 15 
160 

184 ± 18 


146 ± 15 
160 

184 ± 18 


Circinus 


0.011 ± 0.002 


75 ± 20 


75 ±20 


IC 2560 


0.029 ± 0.006 


137 ± 14 


137 ± 14 


Milky Way d 


0.036 ± 0.003 


132.5 


132.5 



As in Hu (2008), NGC3227 may not possess a pseudobulge. 
a The SMBH mass has been updated by Das et al. (2007). 
b The stellar velocity dispersion is taken from Sarzi et al. (2001). 
c The stellar velocity dispersion is from Busarello et al. (1996). 
d The SMBH mass has been updated by Falcke et al. (2009) and 
the stellar velocity dispersions are from Walter et al. (2006). 
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Figure 2. Mass Mbb/Mq versus the mean velocity disper- 
sion cr/rTQ (<to = 200 km s _1 ) for galactic pscudobulges and nine 
GCs together in a log-log plot. The line is the least-square fit 
to the combined data of galactic pseudobulges and nine GCs by 
log(M BH /M ) = 7.3529 + 3.4125 log(o-/<xo) with n = 0.7070. 



nitudes, e.g. Veilleux et al. (1999) reported r c of NGC3079 
to be ~ 3.6 kpc, while Busarello et al. (1996) estimated r c 
of NGC3384 to be ~ 1 kpc and Cavichia et al. (2011) esti- 
mated r c of Milky Way to be ~ 1.4 kpc. 

Fig.[2]shows the correlation between Mbh and the mean 
velocity dispersion a of the host for samples of galactic pseu- 
dobulges and of nine GCs. The joint least-square fit is 

log(M BH /M ) = 7.3529 + 3.4125 log(cr/200 km s" 1 ) 

with n = 0.7070. The n value of A^bh — a power law for nine 
GCs alone is 0.7241. These two n values are close. It seems 
that the Mbh — a power law may extend down to GCs. 



3 RELATION FOR M B n - M GC POWER LAW 

By our model analysis, the total mass Mgc of a GC is 

3n/2 

w,,..= '," . — : 4 I k 3/2 



(3n - 2) 



GC Name 


Af GC (lO 6 M ) 


Major Relevant References 


47Tuc 


1.26 


Pryor & Meylan (1993) 


NGC2808 


1.46 


Servillat et al. (2008) 


lo Cen 


3.1 


Miocchi (2010) 


M80 


1.0 


Pryor & Meylan (1993) 


M62 


0.63 


Pryor & Meylan (1993) 


NGC6388 


2.6 


Lanzoni et al. (2007) 


NGC6397 


0.062 


Heggie & Giersz (2009) 


M15 


0.44 


van den Bosch et al. (2006) 


Gl (M31) 


7.37 ±2.15 


Ma et al. (2009) 
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Figure 3. Masses of central IMBH M B h (in Mq) versus GC 
masses (in Mq) of the nine GCs respectively. The straight line is 
the least-square fit to the data in a log-log display by Mbh/Mq = 
1.25 X 1O~ 4 (M GC /M ) 1 ' 2128 with n = 0.7252. 



By eq a smaller p c corresponds to a larger r c and thus 
a larger A^gc with 2/3 < n < 1. By relation ([7]), Mbh and 
Mgc are related by the power law 



M BH = 



4np c n 
3n - 2 



1/3 



2G 

„2 



(3n-2)/(2-2n) 



M GC 



l/(3-3n) 



In Fig. O we show the IMBH mass A/bh vs the GC mass 
Mgc- By the least-square fit to these data, we obtain 

Mbu/Mq = 1.25 x IO^A/gc/Mq) 1 ' 2128 

with n = 0.7252. The p c value is 4.93M pc~ 3 , grossly con- 
sistent with the data in orders of magnitudes. Note that the 
n value of Mbh — a power law for the nine GCs is 0.7241. 
These two n values are fairly close, indicating that our model 
may consistently explain both the A/bh — a and Mbh — Mgc 
power-laws for the nine GCs. 

In Table 4, we list the central SMBH masses and the 
stellar masses of galactic pseudobulges (Hu 2009) . As noted, 
NGC3227 is not taken as a galaxy having a pseudobulge. 
Moreover, NGC2787 and NGC3384 are two galaxies with 
composite structures consisting of both pseudobulges and 
small inner classical bulges (e.g. Erwin 2008); we do not 
treat them as pseudobulges. We include the Milky Way as 
having a pseudobulge with A/bh = (3.6 ± 0.3) x 1O 6 M 
(e.g. Falcke et al. 2009) and stellar bulge mass M s — (1.3 ± 
0.5) x 1O 1O M (e.g. Dwek et al.1995). In Table 4, M S)B -v is 
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Table 4. A sample of SMBHs in galactic pseudobulgcs 



Galaxies 


logM BH (+,-) 


\ogM a}B _ v 


NGC1068 


7.18 


10.36 


NGC3079 


6.40 (0.30,0.30) 


10.19 


NGC3393 


7.49 (0.03,0.03) 


10.57 


Circinus 


6.04 (0.07,0.09) 


9.57 


IC 2560 


6.46 (0.08,0.10) 


10.26 


Milky Way 


6.56 


10.11 



Here, logAfBH(+,— ) is the logarithm of the mass and lcr error 
of the SMBH and log M s b-v is the logarithm of the host pseu- 
dobulge stellar mass inferred from the ii"-band mass-to-light ratio 
M/L derived from B — V colours. The SMBH mass and the bulge 
mass of NGC 1068 are from Das et al. (2007). The SMBH mass 
of Milky Way is from Falcke et al. (2009). 

the pseudo-bulge stellar mass estimated by K-h&nd mass-to- 
light ratio M/L derived from B — V colour. We do not use 
the pseudobulge stellar mass calculated by K band mass-to- 
light ratio M/L derived from r — i colour, because IC2560, 
Circinus and NGC3393 lack such data. Adopting the least- 
square criterion to the data in Tables 4, we obtain 

Mbh/M q = 4.2825 x lO^M^/Mo) 1 ' 3813 

with n = 0.7587. The p c value is 0.08M Q pc~ 3 . This esti- 
mated p c value appears grossly consistent with the data in 
orders of magnitudes (see Section 2). 



4 CONCLUSIONS AND DISCUSSION 

In summary, for a self-similar quasi-static spherical dy- 
namic evolution of a general polytropic GC (LJ), our model 
analysis shows M B h = £cr 1/(1 ~ n) and Mbh = [(3n - 
2)/(47r Pc n)]( 3 "- 2 '/( 6 "- 6 )[c 2 /(2G)]( 3 ' l - 2 )/( 2 "- 2 'M G c 1/(3 - 3n) 
with 2/3 < n < 1. They agree grossly with current data. 

First, we have the exponent 1/(1 — n) > 3. Secondly, 
n is independent of uncertainties in estimates of r c and p c . 
Thirdly, the tight correlation supports a causal connection 
between the formation and evolution of an IMBH and the 
dynamics of host GC. Finally, while forming an IMBH at 
the GC centre, the spherical general polytropic GC relaxes 
in a self-similar quasi-static phase for a fairly long lapse. 

Index n holds the key in our self-similar quasi-static GC 
model. As p = AK 1/n r - 2/n /(4ttG) and Mbh = Ca 1 ^ 1 -^, 
the smaller the value of n is, the steeper the density profile 
becomes and the smaller the exponent 1/(1 — n) > 3 of the 
Mbh — a relation is. When density profile is less steeper, 
there are more materials around the centre and the accre- 
tion is more effective. Then it would be more effective to 
form IMBHs given that IMBHs are formed by the collapse 
of mostly stars towards the centre. The outcome is that for a 
certain <r, the smaller the mass of an initially formed IMBH 
is, the steeper the density profile is and the smaller the n. 

We conclude that GCs with IMBHs and pseudobulges 
with SMBHs might share qualitatively similar Mbh — cr and 
Mbh — Mgc power-law relations in general. These results 
would bear significance for our theoretical understanding of 
the dynamic evolution of GCs, the formation of IMBHs, the 
connection between GC and pseudobulge formation. 
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